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(57) Abstract: The method of reducing contamination of electrical signals recorded in the presence of repeated interference con- 
tamination comprises obtaining an electrical signal recorded in the presence of a contaminating signal, and detecting a timing signal 
that occurs at a fixed time point during the electrical signal relative to the onset of the contaminating signal. The electrical signal is 
digitized, wherein the digitizing begins with the timing signal. A plurality of digitized electrical signals is analyzed, wherein the elec- 
trical signals are synchronized with respect to the timing signal, to obtain an estimated contaminating signal that is subtracted from 
the digitized electrical signal. This method can be used with electrophysiological signals, such as EEG, ECG, EMG and galvanic 
skin response, and for elimination of noise associated with concurrently used methods such as MRL The method of noise reduction 
is applicable to recordings of other electrical signals, including audio recordings. 
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METHOD AND APPARATUS FOR REDUCING CONTAMINATION OF AN 

ELECTRICAL SIGNAL 

This application claims the benefit of U.S. provisional patent application numbers 
5 60/225,389, filed August 15, 2000, and 60/267,337, filed February 7, 2001, the entire 
contents of each of which are incorporated herein by reference. Throughout this 
application, various publications are referenced. The contents of these references are 
incorporated by reference herein, in order to describe more fully the state of the art 

TP r rH^rr.AT f vrvxp of the invention 

10 This invention relates generally to methods and apparatus for signal processing and data 
collection that are particularly suited for minimizing artifacts and optimizing signal-to- 
noise in simultaneous recording of electroencephalographic (EEG) and Magnetic 
Resonance Imaging (MRI) signals, as well as other environments in which electrical 
signals are subject to repeated interference. The methods of the invention can be applied 
15 to other recordings containing repeated electrical interference, including 

eletromyelographic (BMG), electrocardiographic (ECG) or galvanic skin resistance 

i 

(GSR) signals recording during fMRI and audio recordings or transmissions in the 
presence of 60 Hz noise or electrical transients. 

3ACKGRQUND QF THE INVENTION 

Electroencephalography (EEG) and functional MRI (fMRI) induce mutual artifacts 
when recorded concurrently. Electroencephalography (EEG) has been a key tool for 
study of the brain for decades. However, despite its multiple clinical and research uses, 
such as in epilepsy (Ebersole, 1997), sleep staging (Rechtschaffen & Kales, 1968) and 
psychophysiology, little is yet known about the underlying generators of EEG activity in 
humans. Functional MRI (fMRI) recorded in concert with EEG can provide a method 
for localizing these sources. By using the EEG signal as a reference for fMRI maps, 
concurrent EEG/fMRI opens a new avenue for investigating specific brain function. 
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There remains a need for a system for simultaneous recording of EEG and fMRI, which 
can be used as a tool to localize sources of the EEG. 

Simultaneous recording of EEG and fMRI has proven challenging. Time varying 
5 magnetic (B) fields induce an electromotive force (ejn.f.) in a wire loop perpendicular to 
the B field direction which, by Lenz's Law, is proportional to the cross sectional area of 
the wire loop and to the rate of change of tie perpendicular magnetic field (dB /dt). 
When EEG leads are placed inside the MR scanner, the rapidly changing gradient fields 
and the radio-frequency (RJF) pulses required for MRI may induce voltages that obscure 

10 the EEG signal (Huang-Hellinger, et aL, 1995; Ives, Warach, Schmitt, Edelman, & 
Schomer, 1993). The induced e.m.f. yields currents that can cause heating of the 
electrodes and leads and potentially impart burns to the patient (Lemieux, Allen, 
Franconi, Symms, & Fish, 1997). Motion of the leads themselves within the static field of 
the magnet also induces an e.m.£; even pulsatile motion related to heart beat yields 

15 ballistocardiographic artifact in the EEG that can be of roughly the same magnitude as 
the EEG signals themselves (Ives, Warach, Schmitt, Edelman, & Schomer, 1993; Muri, 
et aL, 1998). Further, introduction of EEG equipment into the scanner potentially can 
disturb the homogeneity of the magnetic field and distort the resulting MR images. 

20 In addition to the large artifacts in the EEG caused by high frequency gradient and RF 
pulses, the high pass filters of most EEG equipment lead to long signal recovery times 
once titte MR acquisition has terminated (Krakow, et aL, 1999). One method used to 
overcome these difficulties in studies of epilepsy has been to monitor the EE G in the 
absence of scanning while the patient is in the magnet and to then trigger functional 

2 5 scanning manually after identification of inter-ictal spikes in the EEG record (Krakow, et 
aL, 1999; Seeck, et aL, 1998; Warach, et aL, 1996). Visual evoked potential has been 
studied using interleaved blocks of ERG and fMRI, where the same stimuli are presented 
in each block (Bonmassar, Aoami, Ives, & Belliveau, 1999). In these methods, EEG and 
fMRI are acquired serially, resulting in protocol limitations and problems with data 

3 o analysis. In the triggered method, relevant changes in the EEG can not be seen during 
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functional scanning. Problems also exist with non-uniform MR image contrast, given that 
Tl saturation typically does not reach equilibrium until 3 to 4 TRs after initiation of the 
scan (depending on TR and effective flip angle). Most often this is handled by ignoring 
images acquired in the first 3-4 TRs, but this then leads to an inherent time delay in the 
5 functional scanning. This could be mitigated to some degree by using schemes that 

correct for the Tl -related intensity differences based on the actual TR (DuBois & Cohen, 
2000; Guimaraes, et al., 1998). In the interleaved method, in addition to the former 
confounds, the EEG and fMRI can not be compared directly. 

io SUMMARY QF THE PSTVENTION 

To overcome the limitations in the prior art described above, and to overcome other 
limitations that will become apparent upon reading and understanding the present 
specification, the invention provides a method of reducing contamination of electrical 
signals recorded in the presence of repeated interference contamination. The method 
15 comprises obtaining an electrical signal, wherein the electrical signal was recorded in the 
presence of a contaminating signal, and detecting a timing signal that occurs at a fixed 
time point during the electrical signal relative to the onset of the contaminating signal 

The method further comprises digitizing the electrical signal, wherein the digitizing 
begins with the timing signal A plurality of digitized electrical signals is then analyzed, 

2 o wherein the electrical signals are synchronized with respect to the timing signal, to obtain 
an estimated contaminating signal. The estimated con taminating signal is subtracted 
from the digitized electrical signal, thereby reducing contamination of the electrical 
signal In a preferred embodiment, the analysis to obtain the es timat e of the 
contaminating signal comprises averaging the electrical signals. In some embodiments, 

25 the analysis to obtain the estimate of the contaminating signal comprises calculating a 
weighted average of the electrical signals. The estimate of the contaminating signal 
be biased towards recent events, for example, by adding the nth electrical signal to a 
scalar multiple, w, of the prior estimate of the contaminating signal and dividing this first 
sum by a second sum obtained by adding the series l+w 2 +w 3 +w 4 ... + ... w 11 . The 
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estimate of the contaminating signal further can be multiplied by a scalar prior to the 
subtracting step. 

The method is particularly suitable for electrical recordings which comprise an 
electrophysiological signal, such as an electroencephalographic recording, an 
5 electromyelographic recording, an electrocardiographic recording or a measure of 
galvanic skin resistance. The method is applicable as well to other types of electrical 
recordings, including audio recordings. In some embodiments, the interference 
comprises interference arising from inductively coupled magnetic fields. The 
interference can also comprise interference arising from alternating current (AQ line 
10 noise. 

One advantageous feature of die method of the invention is that the digitizing can be 
performed at a sampling rate below the Nyquist rate for the contaminating signal. In one 
embodiment, the electrical signal obtained is passed through a low pass filter prior to the 
digitizing, at a frequency of approximately one half of the frequency at which the 
15 electrical signal is sampled. For example, the low pass filter may pass signal frequencies 
of less than about 200 Hz. 

The method can be performed concurrently with Magnetic Resonance Imaging of the 
subject In one embodiment, the electrical signal comprises an electrophysiological signal 
and the contaminating signal comprises gradient activity. Examples of a contaminating 
2 o signal include radio frequency transmitter activity. In a preferred embodiment of the 
method, the di gitizing is performed at a rate of about 200 to about 5000 samples per 
second. The digitizing can be performed at rates below 200 and above 5000 samples per 
second, with representative rates including 100, 250, 500, 1000, 2000, 3000, 4000 and 
6000 samples per second. 

25 The invention additionally provides a method of removing a DC offset from the 

electrical signal by analog subtraction prior to the digitizing. Preferably, the DC offset is 
measured and subtracted from the electrical signal using a difference amplifier. In one 
embodiment, the DC offset is measured by analog to digital conversion, and averaged 
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over a time period long compared to the lowest frequencies of interest in the electrical 
signaL An example of such a long time period is approximately 10 times longer the 
lowest frequencies of interest in the electrical signaL For example, where lowest 
frequencies of interest are approximately 3 Hz, the time period is about 30 seconds. In 
5 one embodiment, the analog subtraction comprises converting the averaged signal to an 
analog voltage and electrically subtracting the averaged signal from the electrical signal 
through differential amplification- In another embodiment, the DC offset is measured in 
an analog integrator having a time constant long compared with lowest frequencies of 
interest in the signal. 

The method of the invention is useful for electrophysiological recordings, such as in an 
electroencephalogram that is recorded concurrently with magnetic resonance image 
acquisition. In a preferred embodiment, the electrophysiological recordings are used to 
inform interpretations of magnetic resonance images. The electrophysiological 
recordings can be used in a statistical analysis of change in intensity of the magnetic 
resonance signaL 

The method can further comprise determining a correlation between change in intensity 
of the magnetic resonance signal and a feature of the electrophysiological recording. The 
correlation can be used to make statistical images, or image maps, that represent an 
association between the electrical signals and the intensity of the magnetic resonance 
signal intensity. In one embodiment, the feature of the electrophysiological recording 
comprises a time course of signal intensity change in defined frequency bands contained 
in the electrophysiological recording. 

The defined frequency bands can be selected to correspond to standard ranges used for 
clinical interpretations of the electroencephalogram. Representative standard ranges are 
selected from the group consisting of from 0 to approximately 4 Hz (the Delta band), 
from approximately 4 to approximately 8 Hz (the Theta band), from approximately 8 to 
approximately 12 Hz (the Alpha band), from approximately 12 to approximately 30 Hz 
(the Beta band), and from approximately 30 Hz and greater (the Gamma band). 
Typically, the frequency bands in this context will not extend beyond 300 Hz. In one 
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embodiment, the method further comprises convolving the time course of the 
electrophysiological signal 'with an estimate of the magnetic resonance hemodynamic 
impulse response function. In this embodiment, the time course of the 
electrophysiological signal is suitably conditioned to more accurately reflect the 
5 anticipated time course of the magnetic resonance signal change. 

The invention additionally provides a method of reducing magnetic interference d uring 
electrophysiological recording from a subject by measuring an electrical potential 
difference between a pair of electrodes, wherein the pair of electrodes communicate with 
a differential amplifier via electrical connections, the method comprising twisting the 

10 electrical connections together, thereby reducing magnetic interference. Also provided is 
a method of reducing magnetic interference during electrophysiological recording from a 
subject by measuring an electrical potential difference between a pair of adjacent 
electrodes, wherein each electrode comprises two leads, the method comprising twisting 
each lead together with a lead of an adjacent electrode, thereby reducing magnetic 

15 interference. In one embodiment, the electrophysiological recording comprises an 

electroencephalogtaphic recording. The method can be performed concurrently with 
Magnetic Resonance Imaging of the subject 

The invention further provides an apparatus for processing digitized electrical signals in 
the presence of a repeated contaminating signal The apparatus comprises a signal 
2 o processor adapted to receive a recording of an electrical signal; a detector adapted to 

detect a timing signal that occurs at a fixed time point during an electrical signal relative 
to the onset of a contaminating signal; a signal accumulator to contain the estimated 
contaminating signal; and a processor adapted to subtract averaged waveforms from an 
electrical signal The signal accumulator can be, for example, a signal averager. 

25 J3RIEF DESCRIPTION OF THg flQURgS 

Figure 1A. Digital photograph of chained bipolar dual-lead dress. Leads of consecutive 
electrodes are twisted together to reduce scanner artifact 
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Figure IB. Schematic diagram showing electrode connectors 2 attached to the head of a 
subject 1 using standard electrode geL Each of the connectors is attached to 2 electrical 
wires (typically constructed of carbon fiber material, which reduces magnetic artifacts). 
The wires from adjacent electrodes are twisted tightly together in pairs 3, where each of 
5 the two wires from a single electrode is twisted together with a different neighbor. The 
electrode pairs are then presented to the input of a differential amplifier 4 where the 
electrical potential difference is amplified to form the electroencephalogram. The inputs 
to the amplifiers are bridged such that the paired leads enclose a complete loop, thereby 
minimising additional differential potentials between amplifiers. The wires are drawn in 
10 thick and thin lines for clarity only. 

Figure 2A. Schematic diagram showing how dual lead electrodes allow each bipolar pair 
3 to be twisted together for their entire length, sending signal directly to local differential 
amplifiers 4. 

Figure 2B. Schematic diagram showing how twisting of leads leaves only small loops at 
15 the head in which e.m.f. can be induced. Current induced in lead twists by motion and 
gradient switching will be self canceling. 

Figure 3. Diagram of an EEG data pathway. KEG signal is fed to a local differential 
amplifier 4, digitized and then sent out of the scanner room 7 via optical fiber for real 
time display and off-line analysis. 

2 o Figure 4. Traces illustrating the ballistocardiogram subtraction algorithm, shown using 
data collected on a normal volunteer. A) A segment of the subject's QRS wave was 
correlated to their EKG, and peak correlation values initiated a trigger (shown by the 
vertical dotted lines). B) Trigger to trigger segments were then averaged, weighting 
segments less and less by temporal displacement from the nth (shaded) segment. Q The 

25 averaged data was subtracted from the raw EEG to yield ballistocardiogram-free EEG. 
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Figure 5. EEG of a phantom recorded inside the MR scanner using twisted and 
untwisted leads, showing recordings both in absence of scanning and during EPI. With 
or without sc anning , EEG recorded using the untwisted leads was significantly noisier. 

Figure 6. Raw EEG of a normal volunteer recorded inside the MR scanner using 
5 untwisted (top) and twisted (bottom) leads. EEG recorded with the untwisted leads was 
significantly noisier. Ballistocardiogram is visible in the twisted lead data, but is reduced 
compared to the untwisted leads. 

Figure 7. EEG recorded on a volunteer during fMRI, before (above) and after (below) 
post-processing to remove the ballistocardiogram. 

Figure 8. Power spectrum of EEG recorded simultaneously with fMRI, in steps equal to 
the TR of 2.5 seconds, showing expected increases in the alpha band (8-12 Hz) when the 
subject's eyes were closed. 

Figure 9. Graph depicting how errors in gradient noise cancellation will occur when the 
sampling and gradient activity are asynchronous. In this Figure, the sampling used to 
create an error estimate (open circles) has drifted by approximately 200 us compared to 
the current sample (closed circles). Subtracting the error estimate actually increases the 
. residual error (by more than 17%), as indicated by diamonds. 

Figure 10. Depiction of timing for the gradient echo EPI pulse sequence used in Figures 
11-13. 

Figure 1 1. Top: Graph representing uncorrected and corrected signals obtained during 
imaging, using 10 kHz sampling. Bottom Enlarged (10X) views of the 25 msec periods 
indicated in dashed lines (uncorrected signal in dotted lines). Note that the artifact 
* suppression varies from cycle to cycle, as a result of phase errors. 

Figure 12. Graph showing simulation of the effects of sampling rate on die efficiency of 
gradient artifact suppression. The trace at the top is actual EEG data recorded at 10 kHz 
during an echo-planar imaging sequence. The three traces below it are the difference 
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between the original signal and a sample lagged at 100 us, 200 us and 400 us as indicated 
(the worst case errors for 10 kHz, 5 kHz and 2500 Hz sampling). The graph below 
shows a detail of tie period during which echo-planar readout occurs. 

Figure 13. Graph similar to that shown in Figure 4, depicting a simulation that shows 
5 the effects of timing errors when sampling at 200 Hz - well below the Nyquist frequency 
for oae gradients. The waveforms are clearly undersampled and therefore appear quite 
different from the prior Figure. The magnitudes of the residual artifacts, however, are 
very similar. The three lower curves show the artifacts that r emain after correction with 
timing errors of 100, 200 and 400 us. 

10 Figure 14. Schematic representation of a low-cost offset nulling differential amplifier 
circuit for use in fMEI. Power supply connections are omitted for clarity. All capacitor 
values are in microfarads. 

Figure 15. MR gradient activity recorded with triggered 200 Hz sampling using a 3s 
repetition time £TR) and nineteen slices. The uncorrected signal is shown at top for a 
15 single TR In the middle is shown the average of 30 TR periods, and at bottom, the 
difference between the uncorrected and averaged signals ("corrected"). 

Figure 16. Human EEG data collected during echo-planar functional imaging. The 
uncorrected data appear at bottom. Above them are corrected records from twenty 
successive TR periods. 

2 0 Figure 17 A. Graph showing energy as a function of frequency and time, derived from 
EEG data acquired during scanning. 

Figure 17B. Graph showing estimated fMDRI activation time course for the EEG data 
appearing in Figure 9A. For clarity, only the Alpha (solid line) and Theta (dashed line) 
bands are shown. 

2 5 Figure 1 8A-C Functional MRI statistical maps of signal change correlated with spectral 
energy at each of five frequency bands (18A left, Delta; 18A right, Theta; 18B left, Alpha; 
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18B tight, Beta; 18C, Gamma), expressed as coefficient of correlation. Note that the 
color scale for the lower frequencies (18A) is different, as the correlations were overall 
lower. All five images were calculated from the same 4:30 (min:sec) acquisition taken 
with the subject at rest with eyes open. 

5 Figure 19. Functional block diagram of MRI-compatible EEG amplifier. 

Figure 20. Alternative differential amplifier circuit for use with fMRL 

I 

Figure 21 . Row chart illustrating method of signal correction by artifact reduction. 
Figure 22. Schematic illustration of elements of an apparatus of the invention. 

DRTATTKD DESCRIPTION 

10 The invention is based on the discovery that contamination of a digitally encoded 
electrical signal can be reduced significantly by making use of a timing signal that is 
associated with the onset of a repeated contamination signal. Such a timing signal can be 
used to align the digitization of repeated contamination signals for determining an 
estimate of the contamination which can then be subtracted from the electrical signal 

15 This method is particularly suited for use with electrophysiological signals, such as EEG, 
ECG, EMG and galvanic skin response (GSR), and for elimination of noise associated 
with concurrently used methods such as MRI. Although the examples described in detail 
herein address the application of the method to recording an EEG in the presence of 
£MRI, those skilled in the art of signal processing will appreciate that the method of 

20 noise reduction is applicable to recordings of other electrical signals, including, for 
example, audio recordings, wherein it is desired to reduce or eliminat e one or more 
sources of repeated contamination. 

All scientific and technical terms used in this application have meanings commonly used 
25 in the art unless otherwise specified As used in this application, die following words or 
phrases have the meanings specified 
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As used herein, "electrical signals synchronized with respect to a timing signal" means 
that data corresponding to electrical signals tecotded over a period of time are aligned so 
that a timing signal that occurs within each of the electrical signals recorded over time is 
superimposed by the alignment When an average is calculated of the electrical signals 
5 superimposed in this manner, the timing signal, as well as any other signal that recurs at a 
fixed time point relative to the riming signal, will be enhanced relative to non-rerarring 
signals. 

As used herein, "twisted" means united by winding, intertwining or coiling. A pair of 
electrical connections between a pair of electrodes and a differential amplifier is 
sufficiently twisted if the enclosed magnetic fields are substantially reduced 

As used herein, "a" or "an" means at least one, unless the context clearly indicates 
otherwise. 

Method? 

The invention provides a method of reducing contamination of electrical signals 
recorded in the presence of repeated interference contamination. The method comprises 
obtaining an electrical signal, wherein the electrical signal was recorded in the presence of 
a contaminating signal, and detecting a timing signal that occurs at a fixed time point 
during the electrical signal relative to the onset of the contaminating signaL The method 
further comprises digitizing the electrical signal, wherein the digitizing begins with the 
timing signaL A plurality of digitized electrical signals is then analyzed, wherein the 
electrical signals are synchronized with respect to the timing signal, to obtain an 
estimated contaminating signaL The estimated contaminating signal is subtracted from 
the digitized electrical signal, thereby reducing contamination of the electrical signaL 

In a preferred embodiment, the analysis to obtain the estimate of the contaminating 
signal comprises averaging the electrical signals. In some embodiments, the analysis to 
obtain the estimate of the contaminating signal comprises calculating a weighted average 
of the electrical signals. The use of weighted averages can serve to achieve an adaptive 
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artifact reduction. The estimate of the contaminating signal can be biased towards recent 
events, for example, by adding the nth electrical signal to a scalar multiple, w, of the prior 
estimate of the contaminating signal and dividing this by the sum of the series 
l+w 2 +w 3 +w 4 ... + ... w a . The estimate of the contaminating signal can be multiplied 
5 by a scalar prior to the subtracting step. 

Figure 19 is a general functional block diagram of representative analog electronics useful 
in the method of the invention. The method perf orms optimally with an adequately 
linear electronic signal. With reference to Figure 19, signal is applied differentially at the 
input terminals 42. Using matched passive components for each lead 44, contaminating 
signal from sources such as radio frequency is attenuated before being differentially 
amplified using standard components 46. The differential amplifier 46 is commonly 
provided with an offset reference input, such that voltages appearing at this terminal ate 
subtracted from the output Using a sample and hold device 52, the amplified DC offset 
potential, derived from the inputs, is first measured, then applied to the differential input 
amplifier 46. To attenuate sources of contamination that are above the highest 
frequencies of interest in the signal, an active low pass filter 48 is supplied. The filtered 
signal is buffered in the output amplifier 50 before being made available to the digitizer 
circuit A switching means 54 is provided to the sample and hold circuit to allow 
detection of the input DC offset at any time desired. 

The circuit modeled in Figure 19 is an example that is especially suited to the problem of 
recording electrophysiological signals during Magnetic Resonance Imaging. The 
functional logic of this diagram is shown in somewhat more detail in Figure 20 which, in 
addition, illustrates a means of transmitting the EEG data out of the MM suite. One 
skilled in the art will see immediately that many different circuit topologies are possible 
that will accomplish essentially the same function for this or other applications. For 
example, in some embodiments it may be desirable to apply a DC offset after the high 
pass filter (as indicated in the schematic of Figure 14), or to avoid this step altogether, if 
significant DC offsets are not present in the signal. Further, in cases where there exists 
adequate headroom for the differential amplifier 46 to remain in its linear range for all 
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expected inputs, it may be desirable to substitute a bandpass filter for 48. The passive 
filter 44 may be eliminated when any signals that reach the input to die input differential 
amplifier 46 are of suitably low amphtude to avoid saturation effects. 

With reference to the diagram shown in Figure 20, EEG signal from the subject 1 is 
5 carried via twisted pair leads 3 to die input of a battery powered head amplifier. The 
inputs include RF attenuation via a series inductance and a parallel capacitor. These 
inputs are mixed in the inputs to shield driver 64, whose output is applied to a concentric 
shield 56 surrounding the twisted leads 3 and connected to the subject 1. The inputs are 
coupled to a differential amplifier 58, whose output is applied to an isolation amplifier, 

10 such as the ISO 122 from Burr-Brown Corporation. This device, and a second ami]** 
device 62 provide electrical isolation to the subject and an added safety factor for EEG 
recording. The output from isolation amp 60 is sampled by an analog to digital 
convertor 66, whose digital output is stored in a latch 68 and converted to an analog 
voltage by digital to analog convertor 70. This output is applied as a DC correction 63 to 

15 differential amplifier 60 after electrical isolation by isolation amplifier 62. 

The output from this head amplifier is presented to a low-pass filter to attenuate signal 
outside of die desired range of the EEG signal This single channel output may be 
multiplexed with the outputs of other similar amplifiers by analog multiplexor 74 clocked 
by hardware timer 78. The output of the multiplexor 74 may be converted to an optical 
2 0 signal by optocoupler 76 and transmitted by .optical fiber, together with the clock signal 
to a second optocoupler 82 that is located outside of the MR shielded room through a 
penetration panel 80. De-multiplexor 82 is used to separate the signals from multiple 
amplifiers and its output may be presented to a differential line driver 88 for transmission 
of the signals over long distances to an analog to digital convertor for later processing. 

2 5 Figure 21 is a flow chart illustrating an application of the method of the invention. The 
raw digitized signal, containing both the desired signal and contaminating artifact are 
shown as (1) in Figure 21. It is strongly preferred, for the invention to work optimally, 
that the signal in (1) be faithfully (linearly) recorded and that the digitization be timed 
with adequate precision to the artifact To produce the corrected signal^ an estimate of 
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the artifact (2) is simply subtracted from the taw signal Optionally, the estimate of the 
artifact may be multiplied by an amplitude constant, m, to account for differences in the 
magnitude of the coupling of the artifact to the signal recording system. 

The artifact estimate can be computed as follows: Each time a new raw sample (1) is 
5 available, it is added to the current estimate of the artifact (2) which has been multiplied 
by a scalar amplitude constant, w (4). This summed signal is then divided by the sum of 
the series 1 + w + w 2 + w 3 4- ... resulting in a new representation of the artifact (5), 
which then replaces the value used for (2). When w is a number less +h» n 1 } the process 
results in a leaky average, where more recent signals have a greater influence on the 
10 estimated artifact (2) than do less recent signals. In this way, the system adapts to slow 
changes in the artifact, if necessary. For the purposes of this disclosure, we have called 
this implementation a 'leaky averager." 

More formally, if R[n] is the n* raw signal collected, A[n] is the estimated artifact for 
collection /?, and S[n] is the artifact-free signal: 



15 



n 

1=0 

and 

S[n]*B[n]-A[nl 



In the leaky averager discussed above, the influence of temporally distant data decreases 
with time. It would be possible to determine w adapuvely by looking at the history of the 
2 0 estimated artifact If the artifact is changing quickly, w should be made smaller (reducing 
the influence of older frames on the current correction). Conversely, w should be large if 
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the artifact is very stable. In the limit, if the artifact cannot change over time, A[n\ should 
be the simple average of all of the samples (w=l). If the algorithm is to be used in real- 
time, A[n] will be the average of all samples until time n. If it is used "off-line" A[»] will 
be the average of all samples both prior to and after time n. 

5 It is possible to determine the value of m adaptively as well, approximating it as the 
amplitude that minimizes the correlation between S[«r] and A[*-/J. 

The method is particularly suitable for electrical recordings which comprise an 
electrophysiological signal, such as an electroencephalographic recording, an 
electtomyelographic recording, an electrocardiographic recording or a measure of 
10 * galvanic skin resistance. The method is applicable as well to other types of electrical 
recordings, including audio recordings. In some embodiments, the interference 
comprises interference arising from inductively coupled magnetic fields. The 
interference can also comprise interference arising from other sources, such as 
alternating current (AC) line noise. 

15 It is particularly advantageous that, in accordance with the invention, the digitizing can 
be performed at a sampling rate below the Nyquist rate for the contaminating signal. In 
one embodiment, the electrical signal obtained is passed through a low pass filter prior to 
the digitizing, at a frequency of approximately one half of the frequency at which the 
electrical signal is sampled. For example, the low pass filter may pass signal frequencies 

20 of less than about 200 Hz. 

The method can be performed concurrently with Magnetic Resonance Imaging of the 
subject In one embodiment, the electrical signal comprises an electrophysiological signal 
and the contaminating signal comprises gradient activity. Examples of a contaminating 
signal include radio frequency transmitter activity. In a preferred embodiment of the 
2 5 method, the digitizing is performed at a rate of about 200 to about 5000 samples per 

second. The digitizing can be performed at rates below 200 and above 5000 samples per 
second, with representative rates including 100, 250, 500, 1000, 2000, 3000, 4000 and 
6000 samples per second. 
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The invention additionally provides a method of removing a DC offset from the 
electrical signal by analog subtraction prior to the digitising. Preferably, the DC offset is 
measured and subtracted from the electrical signal using a difference amplifier. In one 
embodiment, the DC offset is measured by analog to digital conversion, and averaged 
5 over a time period long compared to the lowest frequencies of interest in the electrical 
signal An example of such a long time period is approximately 10 times longer rtmn the 
lowest frequencies of interest in the electrical signal. For example, where lowest 
frequencies of interest are approximately 3 Hz, the time period is about 30 seconds. In 
one embodiment, the analog subtraction comprises converting the averaged signal to an 
10 analog voltage and electrically subtracting the averaged signal from the electrical signal 
through differential amplification. In another embodiment, the DC offset is measured in 
an analog integrator having a time constant long compared with lowest frequencies of 
interest in the signaL 

The method of the invention is useful for electrophysiological recordings, such as an 
15 electroencephalogram that is recorded concurrently with magnetic resonance image 

acquisition. In a preferred embodiment, the electrophysiological recordings are used to 
inform interpretations of magnetic resonance images. The electrophysiological 
recordings can be used in a statistical analysis of change in intensity of the magnetic 
resonance signaL The method can further comprise determining a correlation between 
2 0 change in intensity of the magnetic resonance signal and a feature of the 

electrophysiological recording. The correlation can be used to make statistical images, or 
image maps, that represent an association between the electrical signals and the intensity 
of the magnetic resonance signal intensity. In one embodiment, the feature of the 
electrophysiological recording comprises a time course of signal intensity change in 
25 defined frequency bands contained in the electrophysiological recording. 

The defined frequency bands can be selected to correspond to standard ranges used for 
c lin i ca l interpretations of the electroencephalogram. Representative standard ranges are 
selected from the group consisting of from 0 to approximately 4 Hz (the Delta band), 
from approximately 4 to approximately 8 Hz (the Theta band), from approximately 8 to 
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approximately 12 Hz (the Alpha band), from approximately 12 to approximately 30 Hz 
(the Beta band), and from approximately 30 Hz and greater (the Gamma band). 
Typically, the frequency bands for this application will not extend beyond 300 Hz. In 
one embodiment, the method further comprises convolving the time course of the 
5 electrophysiological signal with an estimate of the magnetic resonance hemodynamic 
impulse response function. In this embodiment, the time course of the 
electrophysiological signal is suitably conditioned to more accurately reflect the 
anticipated time course of the magnetic resonance signal change. 

The invention additionally provides a method of reducing magnetic interference during 
io electrophysiological recording from a subject by measuring an electrical potential 

difference between a pair of electrodes, wherein the pair of electrodes communicate with 
a differential amplifier via electrical connections, the method comprising twisting the 
electrical connections together, thereby reducing magnetic interference. Also provided is 
a method of reducing magnetic interference during electrophysiological recording from a 
15 subject by measuring an electrical potential difference between a pair of adjacent 

electrodes, wherein each electrode comprises two leads, the method comprising twisting 
each lead together with a lead of an adjacent electrode, thereby reducing magnetic 
interference. In one embodiment, the electrophysiological recording comprises an 
electtoencephalographic recording. The method can be performed concurrendy with 
2 0 Magnetic Resonance Imaging of the subject 

Apparatus 

The invention further provides an apparatus for processing digitized electrical signals in 
the presence of a repeated contaminating signal The apparatus comprises a signal 
processor 90 adapted to receive a recording of an electrical signal; a detector 92 adapted 
25 to detect a timing signal that occurs at a fixed time point during an electrical signal 
relative to the onset of a contaminating signal; a signal accumula tor 94 to contain the 
estimated contaminating signal; and a processor 96 adapted to subtract averaged 
waveforms from an electrical signal. Representative variations of the apparatus are 
described in Figures 3, 14, 19, and 20. 
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The signal processor 90 adapted to receive a recording of an electrical signal may be, for 
example, an electronic circuit (IC) consisting of a integrated differential amplifier such as 
a an INAl 14 from Burr Brown Corp., with a DC offeet reference provided through 
measurement by a sample and hold IC such as the LF298 from National Semiconductor 
5 corporation, a low pass active filter and an output buffer all made using standard 

operational amplifier ICs. In this embodiment, the detector can also include a m^na D f 
analog to digital conversion, such as a National Instr ument s NI 6031E installed in a 
personal computet. 

The detector 92 adapted to detect a timing signal that occurs at a fixed time point during 
an electrical signal relative to the onset of a contaminating signal may be, for examp le, 
an optoisolator IC whose output is conditioned using an IC such an LN555 from 
National Semiconductor corporation to produces a TTL compatible trigger signal which 
is then presented to an analog to digital converter such as a National Instruments NI 
603 IE installed in a personal computer. 

The signal accumulator 94 may be implemented, for example, in software in the 'C 1 
programming language as a vector of numbers, or in the National Instruments LabView 
programming language as an array of numbers. Either of these may be executed on a 
personal computer. The signal accumulator 94 can be, for example, a signal averager. 
Other methods, in addition to signal averaging, can be used to generate an estimate of 
the contaminating signal. 

The processor 96 adapted to subtract averaged waveforms from an electrical signal may 
be implemented, for example, as a processing routine in tie f C programming language or 
the National Instruments LabView programming language running on a personal 
computer. Those skilled in the art will appreciate variations on the above examples of 
apparatus elements that will serve the same processing, detecting and accumulating 
functions in accordance with the methods of the invention. 
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Overview of EEG and fM RI in Simultaneous Rp r^Mitig ? nd Mapp ing 

Electroencephalography (EEG) is established firmly as a means to probe changes in 
electrical signals recorded from the scalp that accompany behavioral tasks, and as a 
marker for clinical, cognitive or neural states. Determination of the three-dimensional 
5 localization of the EEG signal is ambiguous because the relationship between the actual 
position of multiple electrical dipoles and the distribution of electrical potentials detected 
at the scalp has no unique solution. Functional Magnetic Resonance Ttnqgttig (fMRI) uses 
modulations in the magnetic resonance signal that depend on variations in blood 
oxygenation to distinguish brain regions whose activity is increased or decreased with 
10 task demands. The following describes a set of solutions to the technical problems in 
simultaneous recording of fMRI and EEG and shows that data from the two methods 
may be combined to create tomographic images indicating brain regions whose activity 
changes as a function of EEG signal intensity in the classically defined spectral bands. 

Sig nifiranr e and Interpretation of the EEG 

Study of the electroencephalogram (EEG) is more than a century old (Caton 1875). The 
phenomenon is highly robust namely that electrical potentials exist at the surface of the 
head that are correlated strongly with ongoing cerebral activity and fluctuate with sleep 
stages (Rechtschaffen and Kales 1968; Buchsbaum, Mendelson et aL 1982; Benca, 
Obermeyer etaL 1992), emotional state (Davidson, Schaffer etaL 1985; Davidson 1988; 
Ekman, Davidson et aL 1990; Lambert and Robertson 1999), attention (Klimesch, 
Doppelmayr etaL 1998; Wrobel 2000), therapeutic drug doses (Loo, Teale etaL 1999; 
Alvarez, Lombardi et aL 2000), traits, such as "aggressiveness" (Fishbein, Herning et aL 
1989) and with circulating levels of a wide variety of drugs of abuse (Cesayirli, Little et aL 
1975; Maykut 1985; Tokunaga, Takeichi et aL 1989; Abraham and Duffy 1991; Mannelli, 
JanM etaL 1993; Bauer, Gross etaL 1997). Despite the considerable history and attention 
to this measurable phenomenon, the origin of the EEG signal and the localization of its 
sources (presumably cerebral), is still not known. The situation is somewhat more 
favorable with evoked potentials (EP ! s, in some contexts known as evoked response 
potentials or ERP's) for which the temporally discrete nature, and the motion across the 
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scalp, of the surface potentials combine to give a better indication of the deep dipole or 
dipoles that generate the signals, but it has been difficult to test directly the relationships 
between scalp EP and brain activity, especially for sources significantly below the cortical 
surface. For example, attempts to localize the generators of the brainstem auditory 
5 evoked response by simultaneous recording with depth electrodes or correlation with 
lesions have been conclusive for only a subset of the waveform components present in 
the signal (Starr and Achor 1978; Chiappa, Gladstone etaL 1979; Achor and Starr 1980b; 
Achor and Starr 1980a; Goldie, Chiappa etaL 1981; Cohen and Britt 1982; Chiappa and 
Young 1985). 

10 When Caton first reported on the resting and task evoked electrical activity of the brain 
. in animals (Caton 1875), he was able to determine that there was a transcortical electrical 
potential that changed during periods of functional activity (sensory s timulati on). Caton 
noted that, "Feeble currents of varying direction. . were generally present between diff erent 
points on the cortical surface. Some years later, Berger noted that scalp potentials could 

15 be recorded in humans with properties similar to those of Caton's cortical potentials 
(Berger 1929) and he soon realized that this electroencephalogram varied according to 
the mental state of the subject (Berger 1930). By 1930 Berger had described what he 
called the alpha rhythm, being relatively high amplitude oscillations in the range between 
8 and 12 Hz that were associated with drowsiness. It is now accepted that alpha activity 

20 is associated with a relaxed, awake state, usually with eyes closed. 

EEG is now a routine and essential test in clinical neurology. It provides diagnostic 
information that cannot be gathered through any other commonly obtainable means. 
Indeed, its indispensability derives from this lack of other routine rlim'rgfl tools to assess a 
broad region of cerebral neurophysiology with high temporal resolution. EE G depicts 
25 moment-to-moment changes in cerebral cortical function, and thus is valuable in any 
clinical context where such information would guide medical decision making. Such 
situations are not limited to electrophysiologic abnormalities, as other pathologic 
processes often affect neuronal function and, thus, impact the EEG. These include 
ischemia, metabolic alterations, mass effect, and infection among others (Markand 1984). 
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Because it represents abnormality in neuronal electrophysiology, epilepsy is a common 
clinic al problem that warrants the use of KEG (Bngel 1984). Despite a central place in 
clinical neurology that has endured over most of the past century, the basis of EEG is 
still understood poorly. The generators of the potentials that sum to create the waves are 
5 understood better (McNamara 1994). However, the interaction between cell populations 
to create electric fields at the scalp has proven difficult to study. 

By performing fMRI during KEG, one may obtain complementary information and 
greater understan ding of EEG and the clinical conditions that EEGs may indicate. It 
appears that solid results that characterize the relationship between scalp potentials and 

10 local brain activity will offer great value in guiding the clinician to isolating specific 

abnormalities. In fact, there are already scattered reports that interictal spike discharges, 
associated strongly with clinic al epilepsy, might be used in combination with functional 
MRI (Warach, Ives etaL 1996; Seeck, Lazeyras etaL 1998; Krakow, Woermann etaL 1999; 
Patel, Blum etaL 1999; Symms, Allen etaL 1999; Schomer, Bonmassar etaL 2000) (or 

15 PET (Henry, Sutherling et aL 1991)) to identify surgically resectable lesions , and 
tomographic localization via MEG has been suggested as a means to guide such 
resections (Stefan; Schneider etaL 1990). The use of fMRI during epileptic seizures has 
also been tested with success by Jackson (Jackson, Connelly et aL 1994). 

20 EEG and Imagin g in Sleep Disorders and Staging 

Based on unit activity, stimulation studies, and lesion studies in non-human mammals, 
brain regions active in non-REM sleep include the anterior hypothalamus, dorsal bulbar 
reticular formation, and nucleus of the solitary tract (Jones 2000). Regions likely 
associated with REM activity, and possibly the generation of wakefulness, include the 

2 5 posterior hypothalamus, ventral mesencephalic pons, basal forehrain, and pontine 

reticular formation. Such localized brain activity should also be visible in humans 
through functional neuroimaging. Both PET and SPECT have been used to examine 
regional activity changes (via cerebral metabolism and blood flow respectively) with sleep 
stage (de termine d electtoencephalogtaphicaHy). These imaging studies have broadly 

3 0 indicated that areas thought to be involved in the active generation of rapid eye 
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movement (REM) sleep are active during this stage of sleep and in non-rapid eye 
movement (NREM) sleep. 

Maquet and colleagues (Maquet, Dive et aL 1990) noted that rates of cerebral glucose 
5 metabolism (rCMRGlc) during NREM are lower, overall, compared to those during 
wakefulness (most notably in the thalamic nuclei), and that REM sleep rCMRGlc is 
comparable to the awake state. Furthermore, the greater the depth of the NREM sleep 
(z.e., the greater the amount of cortical synchrony), the lower the rCMRGlc (Ingvar, 
Baldy-Moulinier et aL 1965; Madsen, Holm et aL 1991; Maquet, Dive et aL 1992). 
10 Especially intriguing is the observation that changes in cerebral blood flow show 
substantial regional heterogeneity. REM-associated increases in rCBF have been 
observed in the pontine tegmentum, thalamus, limbic areas, cortical areas (notably the 
anterior cingulate cortex), and visual association areas, with a decrease in the dorsolateral 
prefrontal cortex, parietal cortex, posterior cingulate cortex, and precuneus (Madsen, 
15 Holm et aL 1991; Madsen, Schmidt et aL 1991; Maquet, Peters et aL 1996; Nofzinger, 
Mintun et aL 1997; Braun, Balkin et aL 1998). Interestingly, rCBF increases have been 
observed in the extrasttiate visual areas, though not the primary visual cortex which may, 
as has been hypothesized by those authors, be indicative of some sort of visual memory 
activation during sleep (Braun, Balkin et aL 1998; Maquet and Phillips 1998; Maquet 
1999). 

As an example of the relevance of imaging during sleep, current studies suggest a role for 
sleep in memory consolidation, based on changes in sleep blood flow as a function of 
daytime activities (Maquet, Laureys et aL 2000). There are also important correlations 
between sleep and a variety of psychiatric disorders (Benca, Obermeyer et al. 1992) that 
might expose sleep physiology as a marker for these problems. But, due to the coarse 
temporal resolution of PET and SPECT (Nofcanger, Mintun et aL 1998) and the 
relatively rapid changes in brain activity during sleep, neither of these imaging methods is 
suited ideally for this purpose. Very recently, there have been reports of the use of fMRI 
in assessment of localized signal changes that take place during sleep (Lovblad, Thomas 
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et aL 1999; Home 2000) that indicate an increase in occipital lobe activity and a decrease 
in the frontal lobes d uring REM, in agreement with the PET findings. However, such 
studies cannot be considered definitive without the incorporation of sleep staging by 
electroencephalographs: means (Lovblad, Thomas et aL 1999). EEG-fMRI provides an 
5 excellent solution. 

The sleep electroencephalogram can be defined by characteristic frequency patterns, and 
brief electrophysiological phenomena such as k-complexes and sleep spindles. 
Furthermore, during REM sleep, there is a descending suppression of muscle tone, 

10 saccadic eye movements, and the loss of the cortical synchrony that is a hallmark of 
NREM sleep. Only with temporal resolving power of fMRI will it be possible to study 
activations associated with these transient events. One application of the invention is to 
assess changes in regional brain activity using fMRI and to correlate such activity to 
classically defined sleep architecture and features. For example, the invention can be used 

15 to seek an un ders tanding the brain activity that underlies the general lack of 

responsiveness to external stimuli, the apparent gating of motor output and dream states. 

Functional MRI (£MR1) is now an established method for the locali2ation of focal areas 
of brain activity, chiefly in humans (Cohen and Bookheimer 1994). Although it is 
assumed that the fMRI signal arises from local changes in blood oxygen content (Ogawa, 

20 Lee etaL 1990a; Ogawa, Lee etaL 1990b; Kwong, BeDiveau etaL 1992; Ogawa, Tank et aL 
1992), this theory has not been subjected to extensive direct testing, and the mechanism 
of coupling between neuroelectrical activity and MRI signal changes is still the subject of 
speculation. Nevertheless, the observed areas of signal increase correlate well with both 
extensive literature on neurophysiology and, more recently, with electro-corticogtaphy 

2 5 derived of humans in surgical settings (Schulder, Maldjian et aL 1 998; Roux, Boulanouar 
etaL 1999; Lurito, Lowe etaL 2000). While there is a growing literature that attempts to 
use the localization power of fMRI to aid in the interpretation of the study of evoked 
responses, there is a striking paucity of reports that attempt to reconcile the findings in 
BEG with functional MRI. 
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The classical or posterior alpha rhythm is found mostly in occipital, parietal, and 
posterior temporal regions (Adrian and Matthews 1934), and first emerges at about 4 
months of age as a 4 Hz oscillation, present with eyes closed and blocked with eyes open. 
The ftequency of this rhythm increases with age, reaching about 8 Hz by age 3, and by 
5 about age 10, it reaches the average adult frequency of 10 Hz (Petersen and Eeg- 

Olofeson 1971). Although recognized since the beginnings of EEG, little is known about 
the functional significance of this alpha rhythm; it reflects essentially a state of relaxed 
wakefulness, and can be used as an indirect measure of brain activation, for increased 
alpha band activity is thought to correspond to decreased activity in underlying cortex 

10 (Shagass 1972). Thus, decreased alpha activity, or stimulus-induced alpha blocking, has 
been termed "event related EEG desynchronization" (Pfurtscheller and Aranibar 1977). 
Davidson and colleagues have shown that alpha asymmetry recorded in anterior regions 
correlates with emotional reactivity, and that these asymmetries appear trait-like in 
waking and in sleep (Petersen and Eeg-Olofsson 1971). Studies in qnimglc have 

15 suggested the thalamus as a possible generator of the alpha rhythm (Petersen and Eeg- 
Olofsson 1971). Lopes da Silva demonstrated significant thalamocortical coherences in 
dogs between lateral geniculate nucleus and pulvinar and the cortex (Lopes da Silva, 
Iierop et aL 1973; Lopes da Silva, Vos et al 1980). Recently, in humans, Iindgren and 
colleagues showed an inverse correlation between EEG alpha power and thalamic 

2 o metabolic rate in normal subjects using PET (Lindgren, Larson et aL 1999). 

Challenges in Combining fMRI and EEG 

Even in the best of circumstances, EEG signals recorded in the clinical environment are 
relatively noisy. The effective input resistances are large, and the signals are small As a 
result, Boltzmann noise limits the ultimate signal to noise ratio. With a typical 5 Mf2 
25 input impedance, the Boltzmann noise will be approximately 1.5 jxV (vm* = <4$kTBR , 
where k is Boltzmann's constant, Tis the temperature, B is the bandwidth and R is the 
equivalent resistance) even over the limited bandwidth of 100 Hz or so used in the EEG. 
Since the scalp potentials are typically only a few uV, the signal to noise ratio (SNR) of 
the EEG seldom exceeds 100:1. Because EEG is often analyzed in the spectral domain 
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over narrow bandwidths, the effective SNR for detection of band-limited signal (for 
example, alpha intensity) may be somewhat higher. 

Beyond the thermal noise limit, other factors further reduce the SNR of EEG. Cardio- 
electric (EKG) and myoelectric (EMG) signals add contamination. Corruption by the 
5 EKG is variable across individuals; though always present, it is usually smaller than the 
EEG, but at times can become comparable in amplitude. The EMG is typically a 
contaminant for only brief periods of head or facial muscle contraction, such as eye 
blink, grimace, etc The methods of the invention actually serve to reduce the EKG 
artifact that would be present whether or not the subject is scanned in the magnet AD 
10 told, these sources of noise reduce the useable dynamic range of the total EEG signal 
dramatically for virtually aE purposes. 

Noise in the MRI Environment 

The noise environment for EEG becomes radically worse when subjects are placed 
inside an MR imaging system. Almost all of the noise sources, however, are coupled 

15 magnetically to the EEG. Several of these are non-biological, including: amplifier noise 
in the shim and field gradient amplifiers; large time-varying magnetic fields induced by 
the field gradients during scanning, and radio-frequency signals generated by the scanner 
for magnetic res onan ce induction. By Faraday's law, the magnitude of the voltage 
induced by these time varying fields is proportional to the first time derhrative of the flux, 

20 and thus to the amplitude of the magnetic field, its first time derivative, and the area 
enclosed by any conducting loop. More specifically: 

em./ = ctyl dt , 

where e.m.f. is the induced electromotive force, and <Pis the magnetic flux. In MR 
imagers these sources cannot be reduced in any practical way. The MR imaging gradients 
25 in state-of-the-art ima ging instruments slew at an extremely high rate; the field gradients 
on a typical sc anne r operate at 80 T/ sec and are thus major sources of noise; the newest 
generation of MR instruments, with local head gradients, will slew two to three times 
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faster. The RF pukes, though only about 50 milli Gauss, have fundamental frequencies, at 
the 3 Tesla operating field, of 128 MHz. Slewing at 4000 T/s, they too ate large sources 
of noise. 

Physiological noise sources are present also. Even small motions of the subject are 
5 coupled to the EEG as the leads move within the large static magnetic field. Not only 
minor fidgeting, but also the motion of the whole body with each heartbeat (the 
ballistocardiogram), produces signal in the y.V to mV range. 

DC Offsets and Transient Recovery 

The scalp electrical potentials used in the EEG contain both time varying and static (DC) 
10 components. Often the DC signal is much larger than the EEG, but it is seldom of 

interest for clinical diagnostic purposes (note reference to true DC offsets, not to slowly 
varying potentials), as it contains essentially no information. However, it does cause 
trouble for the EEG in several ways. Typically, the DC offset increases the dynamic 
range needed to digitize the EEG signal For example, the EEG signal maybe only a few 
15 juV, while potentials of a few mV may exist between electrodes, or as a result of the 

chemical electrode potential. The signal digitization depth will be reduced by the ratio of 
the EEG to the DC potential. For example, assuming a DC offset of 10 mV and an 
EEG signal of 10 uV, the 4096 different levels representable by a 12 bit analog to digital 
converter (ADC) will be reduced to only 4 levels for the EEG. Clearly this loss is 
2 0 unacceptable, as quantization noise will dominate the signal. 

For these reasons, conventional EEG amplifiers are equipped with AC-coupled (high 
pass) inputs, usually a capacitor separating the first stage amplifiers from the input to the 
ADC. The inputs will usually have time constants of several seconds, allowing 
frequencies of 1 Hz or so to pass unattenuated. One consequence of this AC coupling is 
2 5 that it creates a time constant for signal recovery if the input saturates. Because these 
filters must pass very low frequencies, the recovery time for the analog signal to come 
back to the center of its nominal range can be quite long. 
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In fMRI, the recovery time associated with AC coupling is a substantial problem, as the 
gtadient-induced artifacts (tens of mV) can be large enough to saturate the input stages, 
pinning them to the positive or negative supply rails for several msec. When the 
gradients cease, the settling time of the high pass filter greatly outlasts the gradient event 
5 A recent paper studying EEG-fMEI combinations, reported that, <c The EEG could not 
be interpreted during the artifact caused by the excitation pulse, but the recording 
becomes readable in less than 1 second (approximately 100 msec) after completion of the 
BURST" (Hennig and Hodapp 1993; Lovblad, Thomas et aL 1999). This problem can be 
mitigated by using the methods of the invention, including an artifact-reducing electrode 
10 configuration and an input amplifier (Grass-Telefactor) with enough headroom to stay 
out of saturation prior to the high pass stage. As disclosed herein, the invention further 
provides a more complete solution that avoids the high pass filter completely. 

Gradient Noise 

The magnetically-induced gradient noise is of very large amplitude (milliV olts in a typical 
15 s cann er) as compared to the EEG, especially in the context of echo-p lanar imaging . One 
group has implemented a correction scheme for gradient artifacts that is similar to a 
scheme that gtoup developed for ballistocardiogram removal (Allen, Josephs et aL 2000). 
Because the fundamental frequencies of the gradient activity are much higher, they 
developed special recording hardware that allowed them to use a much higher digital 
2 o sampling rate of 5 kHz, which they selected as being rapid compared to the nominal 
Nyquist frequency for the gradient waveshapes. Unfortunately this is not sufficient, as 
sampling at the Nyquist rate guarantees only against aliasing of the higher frequencies 
into the pass band, but does not effectively remove the artifact 

For example, assume that there is an undesired signal consisting of a sinusoid at 1000 Hz 
2 5 contarruriating the EEG, which is sampled digitally at 5 kHz. Because the gradient and 
digitization activity ate clocked independently, and are asynchronous, the phase at which 
the artifact is sampled can vary by as much as 2%/S (72° or 200 us at this frequency). 
Over extended sarnpling periods (typically five minutes or more in an imaging 
experiment), it is likely that the relative timing of the scanner and sampling clocks will 
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differ to this degree. Figure 9 shows this effect In Ads simulation, the artifact is assumed 
to have been sampled at the points indicated by closed ci r cles. At a later time, the 
sampling has drift ed with respect to the gradients by 200 us (dashed line, open circles). 
The difference signal is the residual artifact, in this case just over 17% greater than the 
5 uncorrected signal To mitigate this problem, the Allen group has adopted a sophisticated 
interpolation scheme that is successful in tninitni^ing the residual contamination. The 
present invention provides a much more effective approach based on an alternative 
formulation of the digital sampling problem. 

The methods disclosed herein have made it possible to effectively eliminate 
10 contamination of the EEG signal by the most severe sources of noise present during 
MRI s canning in general, and in functional MRI in particular. In one embodiment 
described herein, the method has been used in the construction of tomographic maps of 
brain activity corresponding to the energies in spectrally-defined components of the 
EEG. 

is Theory 

Digitization 

The artifact from gradient activity is large and contains substantial energy at high 
frequencies. Figure 10 shows the timing of a typical echo-planar imaging sequence, as 
used in typical functional studies. The lines for Select, Phase and Readout indicate the 
20 amplitudes of the three orthogonal magnetic field gradients used for imaging. The fourth 
line indicates the timing of the radio frequency channel (only the amplitude envelope is 
shown for the RF, as the carrier frequency of 128 MHz is not visible at this resolution.) 
Immediately apparent is the very large high frequency (1400 Hz) oscillation of the 
readout gradient (shown at half the vertical magnification of the other gradients). 

25 Figure 11 shows raw signal, recorded from an EEG system at a sampling rate of 10 kHz, 
following analog low-pass filtering at 100 Hz. The insets on the bottom of this figure 
show expanded representations (10X) of the indicated regions of the signals. When EEG 
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data are acquired during s cans , the MRI field gradients induce voltages much higher than 
the cortical signal. Comparing this and the previous figure reveals that the magnitude of 
the high frequency components is reduced dramatically by the low pass filter, and that 
the residual artifacts outlast the gradients themselves. 

Though the low pass filter provides at least 100-fold reduction in the 1400 Hz 
oscillations, it does not remove the trail stents as these gradients turn on and off. These 
contain energy at very low frequencies, as welL The AC coupled input stage, in this 
traditional amplifier, is responsible for the extended ring-down of the artifacts (the 
saturation recovery alluded to previously), although it is much better than the 0.1 to 1 s 
ring down reported by others (Lovblad, Thomas etaL 1999), presumably due in large part 
to the attenuation provided by the differential recording apparatus, which helps to 
prevent the amplifiers from going into saturation. 

When applying cyclic averaging techniques to this data set (Allen, Josephs et aL 2000) as 
shown in Figure 11, they are reasonably effective in attenuating the effects of the low 
frequency components, they are largely ^effective at removing the high frequency 
content This, as described above, results from the asynchronous sampling. The residual 
(worst case) error from sampling too slowly can be predicted from the maximum phase 
shift, <p 

e = cos(2tc/? + <p) - C0S(27f/?) 
= cos(27r# )(cos sin(27ty?) sin <p 

20 The maximum phase shift, <p, that can occur at a given sampling rate is equal to 2nfo/ 6, 
where ft is the sampling frequency, and is the frequency of the EPI readout- 
Comparing the residual artifact during the two scan periods (two expanded frames at 
bottom of Figure 1 1) reveals that the cancellation efficiency is unstable as a consequence 
of the asynchronous riming of the gradient activity and sampling device, which causes 

25 the sampling offset, cp, to drift over time. 
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As the sampling rate is increased, the cancellation will become more accurate. Using the 
approximation that for small a, siaa « a and cosa ~ 1, one can see that in this regime e is 
approximately proportional to cp. It follows that if the artifact must be suppressed by 100 
fold, the signal must be digitized at approximately 100*2*7t times tie highest frequency . 
5 of interest, in this case (the 1400 Hz readout) about 880 kHz/channel, which is 

impractical for reasons of both cost and overall data handling. In any case the acceptable 
sampling error is predicted readily by this formula given the low pass filter characteristics 
and the desired final signal to noise. 

The data above were sampled at 10 kHz, well above the 2800 Hz Nyquist criterion for 

10 the high frequency components of the signal. Figure 12 shows the efficiency of artifact 

subtraction based simply on rapid sampling. On the top is shown the raw artifact Below 

it is shown the residual artifact that remains after subtraction if the riming of the 

sampling and the scanner have drifted from synchrony by 100, 200 and 400 usee, 

corresponding to the worst case errors for sampling at 10 kHz, 5 kHz and the 

15 approximate Nyquist rate of 2500 Hz, respectively. It is immediately apparent that the 

residual artifact, after subtraction, is large even with the smallest timing offset. In the 

graph at the bottom of Figure 12, which shows in greater detail only the echo-p l a na r 

readout segment of that data, one can see that the simple subtraction actually increases 

the magnitude of the artifact, as predicted in the equation above (and in Figure 9). 
» 

2 0 The efficiency of the subtraction of the gradient artifacts is effectively independent of 

sampling rate, and a repeated single sample, properly timed, can be used to correct for the 
artifact at that time point completely. This general finding is shown in Figure 13 which 
shows the same effects of timing shifts at a low sampling rate of 200 Hz, well below the 
Nyquist rate (the gradient activity is aliased into the digitized EEG signal). As in the prior 

2 5 example, the magnitude of the residual artifacts increases as the sampling is delayed with 
respect to the scan timing. Its magnitude is no worse than that seen with mote rapid 
sampling. 

It is clear that if the sampling is timed precisely to the gradient activity, the residual errors 
will be eliminated much more effectively. Perhaps less intuitive is the tact that this does 
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not requite Nyquist rate sampling for the artifact frequencies. Time (phase) shift and 
frequency can be seen as "duals": Perfect sampling needed for artifact subtraction could 
be achieved by sampling at infinitely high frequency with arbitrary timing or with precise 
synchronization at an arbitrarily low sampling rate. To yield the same 100:1 suppression 
5 that requires 880 kHz sampling, the sample timing would require an accuracy of 1/880 
kHz, or 1.14 jusec. Ideally, the residual scanner artifact should be small compared to the 
thermal noise of the EEG signal With analog filtering and proper recording technique 
' (outlined in the examples below), the scanner artifact is about ten times the amplitude of 
the EEG. Ass uming an EEG signal to noise ratio of 100:1, a thousand-fold suppression 
10 is needed in the digital processing, achievable with either 8.8 MHz/ channel sampling 

(which could come only at tremendous expense in the digitization hardware), or with an 
easily achieved 1 14 nsec timing accuracy 

EXAMPLES 

The following examples are presented to illustrate the present invention and to assist one 
15 of or dinar y skill in making and using the same. The examples are not intended in any 
way to otherwise limit the scope of the invention. 

Exat nplft 1 : Arrpii^ntr simultaneous EEG and functional MKT 
Methods 

EEG Device and Lead Placement 

2 0 The EEG device incorporates numerous hardware modifications to reduce artifact in 
concurrent EEG/fMRI, and was provided by Tele&ctor Corporation (W. 
Conshohocken, PA). Signal is detected from the scalp using silver chloride plated plastic 
cup electrodes connected to a compact magnet-compatible local amplifier (headbox) via 
10 foot carbon fiber leads with a resistance of 1 k£2/foot This design minimized both 

2 5 artifact in the MR images and the induction of RF current loops in the lead wires. 

A lead configuration was devised that minimized unwanted current induction by 
recording EEG in a hard-wired montage using special dual lead electrodes. The lead 
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wires from each electtode pair were twisted together over their entire length, fo rming a 
chained bipolar montage for each hemisphere (Jp2-j8 } j8-t4, t4-t6, t6-o2 } o2-p4,p4-(4, <4-fl, 
j*-$2;jp1-p,f7-t3 s t3-t5, t5-o1, o1-p3,p3-c3, c3-fi,p-ft>1). The dual leads allowed each 
differential pair to be twisted together (Figure 1A-B and Figure 2A). The resulting 
5 configuration leaves only small loops at the head in which current can be induced. As 
shown in Figure 2B, currents induced by motion and gradient switching should be self 
canceling, for the current induced in consecutive twists will flow in opposing directions 
in each lead wire. 

The magnet-compatible headbox contains 32 separate channel inputs, each with a 
10 differential amplifier coupled to an RC filter having a time constant of 0.25 msec. 

Sixteen channels were used to record EEG, and two additional charmplg were used fot 
electrocardiogram (EKG). EKG was acquired using a pair of twisted single lead 
electrodes placed above and below the heart on the subject's back. This placement 
minimized lead motion, and thus electrical artifact, due to breathing. A scan trigger 
channel was also used to receive a pulse from the scanner every TR to aid in post- 
processing of the data. The signal was filtered in all channels with a band pass of 0.5-70 
Hz to further attenuate high frequency noise. A single lead was connected at to 
headbox ground as an added patient safety measure, but was not used as a montage 
reference. 

The signal was fed to a battery powered isocoder containing an A/D converter where it 
was sampled at 200 Hz., and the digitized signal was c arrie d out of the shielded magnet 
room via optical fiber to maintain the scanner's electromagnetic isolation. After 
translation to TTL, the data were routed to a Telefactor Digital EEG (D/EEG) (486 
computer). The EEG data could then be viewed in real time and sent off-line via a 
lObaseT Ethernet connection to a post-processing and viewing station for further artifact 
attenuation. 
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Scan Protocol 

All subjects signed a consent form approved by the UCLA Human Subject Protection 
Committee prior to MR scanning, which was performed in a General Electric 
(Waukesha, WI) 3T scanner modified for Echo Planar Imaging (EPI) by Advanced NMR 
5 Systems. No visual or auditory stimulation was provided to the subjects during 
functional scanning. 

Scout scans of the entire brain were first acquired to localize slice planes parallel to the 
AC-PC line through the occipital cortex. To acquire EEG during functional scanning, 
the scan protocol was then specified to allow windows of readable EEG between 

10 gradient bursts. An EPI sequence was used with TR = 4000 ms, echo time (TE) = 45 
ms, 64 x 64 matrix, 20 cm x 20 cm field of view (FOV), 4 mm slice thickness, and 1 mm 
gap to collect 6 slices spaced evenly over the TR. period, leaving a 580 ms window of 
readable EEG between each 90 ms period of gradient induced noise (an 87% duty cycle). 
An EPI sequence (TR = 6000, TE = 54, 128 x 128 matrix, 20 cm x 20 cm FOV) was 

15 then acquired coplanar with the functional scans for use as an anatomical reference. 

Artifact Redu ction Post Processing 

Post-processing and viewing was performed on a Dell Inspiron 3000 Pentium PC. After 
importing data from the D/EEG over the wire, the EEG data was viewed with 
Telefactor Twin software, and processed further using home-built software described 
2 0 below to remove remaining artifacts. The residual artifact included both noise from 

e.m.f. induced by the magnetic field gradients (which appeared in the EEG when a slice 
is acquired) and ballistocardiogram. The latter occurred in a fairly regular pattern just 
lagging the Q wave of the EKG, but its morphology and amplitude differed in each 

EEG rh pnnftl. 

25 To suppress gradient noise and ballistocardiogram in the EEG record, the gradient noise 
was first removed by blanking the EEG and EKG data for the duration of each MR slice 
acquisition. Following a trigger every TR, the 90 ms data segments containing scan 
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artifact were replaced with zeros. Therefore, large deflections caused by the MR 
gradients did not corrupt die EEG data during the averaging and subtraction in further 
artifact removal. 

Next, using a method similar to that of Allen et aL (Allen, Polizzi, Krakow, Fish, & 
Lemieux, 1998), sections of EEG were averaged together following a cardiac trigger to 
yield the ballistocardiographic artifact To identity the initiation of each cardiac pulse, a 
single artifact-free QRS wave segment of the subject's EKG, recorded inside the scanner 
when no scanning was taking place, was used as a reference (see Figure 4A). This 
reference wave segment was compared to portions of the EKG data of the same number 
of data points as the reference segment that were shifted by one point at a time, and 
calculated a correlation coefficient (CQ for each data portion. When the CC exceeded 
an empirically selected value (typically 0.7), the peak CC following this threshold crossing 
was identified to trigger the initiation of ballistocardiogram averaging and subtraction. 

Figure 4B illustrates the averaging and subtraction algorithm performed on the data in 
each EEG channel Every trigger-to-trigger section of raw EEG data (An) was averaged 
with all preceding sections. Because EEG and EKG should be uncorrected, this 
method averaged out the EEG signal and left only ballistocardiogram <Bn). Data 
sections were weighted inversely with their temporal displacement from the current 
sample to compensate for slow changes in the ballistocardiographic artifact, calculating 
the ballistocardiogram in each trigger to trigger section using the weighted average 

1-0 

with a weighting factor w = 0.9. Thus, the earlier cycles formed an exponentially 
decreasing contribution with a time constant of roughly 10 sections. The averaged wave, 
B„, calculated separately for each channel, was subtracted from that channel's raw EEG, 
2 5 An, to yield artifact free corrected EEG. 
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Variations in section length due to changes in heart rate were accounted for by averaging 
sections point-by-point. The first data point in each section was averaged with die first 
data point in previous sections, the second with the second, etc. Data points at the end 
of longer sections were averaged with corresponding points in other long sections, and 
s the weighting factors of each point were adjusted accordingly. To avoid subtracting data 
with too few points averaged, the subtraction was performed only if three or more 
points were used in the calculation; otherwise, the raw data remained in the final record. 

rii Q «»rt<>rittali on of Noise Reduction 

Kfnise Reduction due to l**d Dress: Phantom Study 

10 To characterize noise reduction due to the twisted dual-lead dress, scanning experiments 
were performed using a biological phantom - a 9 pound head-sked grocery store 
roasting chicken - and compared twisted vs. untwisted lead arrangements. The leads 
were placed at distances corresponding approximately to standard international 10-20 
positions on the phantom. The eight electrode chained twisted montage detailed above 

is was placed on the left hemisphere, and on the right was placed a matching dual lead non- 
twisted montage. A scout scan was acquired to position the functional dices to cover the 
phantom Functional EPI scans wete then performed as described previously. 

The EEG data were analyzed to quantify noise reduction due to the lead dress, 
calculating the loss, in dB, using 30 second data segments as 



20 



dB ■ 20 Iol . 



where averages were taken over all twisted and all untwisted channels. Gradient noise 
reduction was calculated by subtracting the square root of the sum of the squared 
voltages without scanning from that during scanning, and then calculated attenuation as 
above. 
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Tsfoke Reduction due to Lead Dress: Human Study 

To characterize twisted lead noise reduction further, the above study was repeated on a 
25-year-old normal male volunteer. Again, the chained twisted lead set was placed on 
the subject's left hemisphere, and the untwisted set on the subject's right. The subject's 
5 EEG was recorded inside the MR scanner, both with and without scanning. 

gpecttal Analysis 

By restricting the timing of the EPI acquisitions to fell outside of the frequency band of 
interest for EEG (e.g., to study alpha activity, the sc anning rate must be less than four 
images/second), it was possible to retain useful EEG spectral information. To illustrate 
10 this, the 25-year-old normal volunteer was scanned during three different tasks known to 
moderate alpha power - a basic eyes open/eyes closed task, a math task, and a 
visualization task. All studies were performed during functional MRI as described above. 
A baseline eyes open scan was acquired, and then the three tasks were run as follows. 

In the eyes closed task, the subject was given verbal cues during the scan to keep his eyes 
15 open for the first two minutes, closed the next two minutes, then open for die last 

minute. For the next two tasks, the subject was told to keep his eyes closed. First, the 
subject was instructed to count backwards by sevens from the four digit number given to 
him two minutes into the scan. To verify that the task was performed, the subject was 
asked for the number he ended on after the three minutes of counting. Second, the 
2 o subject was instructed to visualize eating his favorite meal, again two minutes into the 
scan. 

After post-processing die EEG to remove artifact, the alpha power in each TR was 
calculated using software developed in house. With this software, the EEG power in 
user-defined bands with each TR was found using the Fast Fourier Transform. Spectral 
2 5 power in the alpha band was used as a reference function to calculate fMRI signal maps. 
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Results 

Kfnise Reduction due to Lead Dress; Phantom Study 

EEG recorded on the phantom using the untwisted leads was substantially noisier than 
that recorded using their twisted counterparts. Figure 5 shows the recorded EEG in both 
5 twisted and untwisted lead channels. When not scanning, the twisted leads reduced 
random noise power by an average of 5.4 dB. During scanning, gradients caused large 
artifacts in both the twisted and untwisted lead sets, but this noise power was reduced by 
an average of 6.3 dB in the twisted leads. 

Twisted vs. Untwisted L eads on a Volunteer 

10 Figure 6 shows EEG data recorded on a normal volunteer inside die scanner when no 

s canning was t^Tnng place. Twisting the leads (shown here on the left hemisphere vs. the 
untwisted on the right hemisphere) reduced noise by an average of 7.5 dB across all 
channels. 

^fnisa Reduction Due to Artifact Post-Processing 

15 Figure 7 shows EEG recorded on a volunteer during functional MRI before and after 
post-processing. The post-processing removed significant gradient and RF artifact, as 
well as ballistocardiogram. 

Spectral Data 

The data shown in Figure 8 was recorded on a normal volunteer during an eyes open, 
2 0 eyes closed paradigm using simultaneous EEG/fMRL Here, four slices were acquired 
with a 2.5 second TR. When the subject's eyes are dosed, power in the alpha band, 
between 8 and 12 Hz, increases significantly. And as is expected, this alpha signal is not 
present when the subject's eyes are open. 
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Discussion 

Using a combination of analog pre-preprocessing and digital post-processing, this 
method allows one to record clean EEG during functional MR scanning. While the 
EEG is obscured during gradient bursts, it recovers quickly after each slice acquisition. 
5 With this method, then, a trade off must be made between brain coverage in functional 
scanning and the fraction of useable EEG in the data record. Consideration must also 
be given to the timing of slice acquisition so that it does not overlap with desired spectral 
frequencies. 

The invention provides a potentially powerful tool for localizing the sources of various 
10 EEG waveforms. Because the EEG is acquired simultaneously with fMRI instead of 
serially, it can be used as a direct source for the fMRI reference function. In this way, 
activation maps could be made of any relevant changes in the EEG: spike and slow wave 
patterns in epilepsy, spectral changes, or even event related potentials. 

Rumple 2; Method for removal of artifacts in simultaneous EEG and fMRI 

15 The strategy outlined below is applicable to any method involving digital subtractive 
noise cancellation, including EEG, EEG/fMRI, as well as any environment in which 
subtracted averages of artifacts are used to reduce noise. 

Methods 

KRG recording and electrode placement 

2 0 The basic approach of differential recording has been described previously (Goldman, 
Cohen et aL 2000; Goldman, Stern et al 2000). Paired silver electrodes are placed on the 
scalp surface attached with conductive electrode gel and check to ensure that the 
nominal impedance is less than 5 Hi for each electrode. The electrodes are themselves 
attached to two carbon fiber conductors having a distributed resistance of approximately 

25 3 kfi/m. The wires are dressed in pairs with the leads from adjacent electrodes twisted 
together tightly to minimize electromagnetic pickup. This configuration offers about 6 
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dB attenuation of artifacts from gradient and ballistocardiographic noise sources 
(Goldman, Stern et al 2000). 

Amplification 

Further improvements in the overall performance of in-magnet EEG can be achieved 
5 through the use of better analog electronics. Although the ballistocardiogram contains 
substantial energy in the range of most interest to clinical EEG (from 1 to 50 Hz), the 
other ynnin sources — gradient noise and RF transmit noise — have much h igher 
fundamental frequencies. In a typical scanner, for example, the overwhelming majority of 
the gradient-related noise is at a fixed frequency of 1400 Hz, with significant energy 
10 distributions down to 100 Hz or less. The radio frequency eneigy, of course, is at 
radically higher frequencies well outside of the interesting pass band for EEG. 

Disclosed herein is a very simple circuit that corrects for much of the analog portions of 
the artifact (Figure 14). An initial gain stage features a single pole filter to attenuate the 
large RF signals prior to sam pling and provides enough gain to bring the EEG signal into 
15 the mV range without the artifacts causing saturation. The next stage offers 30 

dB /octave attenuation at a comer frequency of 200 Hz (so that gamma range EEG is 
readily passed.) 

The AC coupling problem is handled in the final amplification stage, which is arranged 
to include a resettable offset-nulling circuit that stores any DC offset across a low-leakage 

2 0 capacitor. The offset nulling switch can be, e.g., a mechanical switch, or, so that the 

nulling can be performed tinder digital control as needed, presumably when the software 
detects that the signal is close to digital saturation, a CMOS switch. One can reduce this 
design to use with a pc and assemble a single channel of this system. Under lab test 
conditions, the circuit is able to hold the DC ofbet to within 1% for 10 minutes at a 

25 time, easily meeting the system requirements for MR scanning. 

The AC coupling problem is handled in the final amplification stage, which is arranged 
to include a resettable offset-nulling circuit that stores any DC offset across a low-leakage 
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capacitor. The offset nulling switch can be, e.g., a mechanical switch, or, so that the 
nulling can be performed under digital control as needed, presumably when tlie software 
detects that the signal is close to digital saturation, a CMOS switch. This design has been 
reduced to printed circuitry and a single channel of this system has been assembled One 
5 can reduce this design to use with a personal computet and assemble a single channel of 
this system. Under lab test conditions, the circuit is able to hold the DC offset to within 
1% for 10 minutes at a time, easily meeting the system requirements for MR scanning. 

Imaging 

Presented herein are examples of activation mapping data from two subjects (one each 
10 for the raw EEG data and for the EEG energy maps). Both subjects were without 

neurological or radiological abnormalities as assessed by a brief neurological inventory 
based on a form developed by the National Institutes of Health NIH and a neurological 
inventory performed by a board-certified neurologist Beyond lying in the magnet during 
scanning with eyes closed, the subjects performed no explicit cognitive task. 

15 All sc anning was performed on a General Electric (Waukesha, WT) 3.0 Tesla Signa® 

scanner modified by Advanced NMR Systems (Wilmington, MA) for high performance 
echo-planar imaging (Brady, Cohen et aL' 1991; Cohen, Kelley et aL 1996). For testing of 
artifact rejection, a 19 slice echo-planar data set was collected with a TR of 3 seconds, 
TE =45 ms, 4 mm slice thickness and 3.125 trim in-plane resolution (64 x 64 scan matrix 

20 and 20 cm FOV), to achieve appropriate weighting for "BOLD" contrast effects 

(Ogawa, Lee et aL 1990a). For the mapping data, after image-based shimming (Reese, 
Davis et aL 1995) and collection of a scout scan, imaging was performed using a gradient 
echo EPI scan as described above. These data, however, were acquired with a longer TR 
of 4 seconds and only four slice planes. The pulse sequence was modified to include a 5 

25 us trigger pulse at the beginning of each TR period, the leading edge of which was used 
for synchronization of the EEG sampling acquisitions. 
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Pigitization of the EEG 

EEG data ate acquired using a PCI-1200 (National Instruments, Houston, TX) on a pc- 
compatible microcomputer. Using LabView (National Instruments, Houston, TX), 
sampling software was developed that responds to the leading edge of the scanner trigger 
5 by acquiring a fixed number of samples at a user-specified rate. Specifi c ally, with a TR of 
3 seconds, one can acquire 599 samples at a rate of 200 Hz following each trigger. These 
data are immediately flushed to a file in a background process; to accommodate this 
process one sample point is dropped with each TR. 

Averaging and Artifact Removal 

10 A fame of EEG data is collected from each channel with each scan TR, and the frames 
are averaged separately for each channel to create an accurate representation of the 
s canne r artifact One can then remove the artifact by simply subtracting the averaged 
signal from its respective channel. The next step is to manually inspect the EEG data for 
characteristic artifacts, such as eye blink and facial muscle movement that are recognized 

15 easily by their morphology within and across c h annel s . 

Ballistocardiogram Suppression 

Details of the ballistocardiogram suppression procedures have been published previously 
by Goldman et aL, (Goldman, Stem et ai 2000) and will only be summarized here briefly. 
Because the cardiac-induced motion is nearly repeated with each heartbeat, the resulting 

2 0 artifact is essentially the same and is superimposed onto the desired EEG. To remove 
this artifact one can detect the heartbeat using the electrocardiographic signal and 
calculate the average ballistocardiogram over many heartbeats. One can then subtract this 
average from the EEG signal. Of course, die morphology of the ballistocardiogram 
differs in each of die EEG channels as the motion of each lead differs slightly. Therefore 

25 this averaging and subtraction process is preferably performed separately for each of the 
EEG lead pairs. The method is conceptually like that described by Allen et aL, (Allen, 
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Polizzi et ai 1998) but differs in the details by which slow variations in the 
ballistocardiogram ate accommodated 

fcnage Processing 

Using software developed in-house, standard FFT methods were used (Press, Vetterling 
5 etaL 1992) to determine the power separately in five different spectral bands in the EEG 
signal for each TR period Following convention, data were combined from 0.5 to 4 Hz 
("delta"), 4 to 8 Hz C^eta' 3 ), 8-12 Hz ("Alpha"), 12 to 30 Hz ^Beta") and from 30 to 
70 Hz ("Gamma"). Using this model of spectral energy as a function of time, one can 
then estimate a prediction of the BOLD signal changes by convolution of the spectral 
10 data with an a priori model of the brain hemodynamic response function (Cohen 1997). 
The convolution introduces a lag in the time course that is thought to represent a 
reasonable estimate of the hemodynamic latencies and, in addition, acts as a low pass 
filter, which tends to reduce somewhat any noise in the EEG data (see Figure 17). 

One can then use scanSTAT (available via the internet at the URL for the UCLA Brain 
15 Mapping Center, http:/ /www.brainmapping.org) to first spatially filter the images and 
then to form statistical maps indicating the correlation between the EEG spectral power 
in each band and the local fluctuations in MR signal intensity. The convention is to 
indicate regions of increasing positive correlation in colors from red to yellow and areas 
of increasing negative (anti) correlation in colors from blue to cyan. In particular, regions 
20 of high negative correlation ate interpreted as indicating decreased blood flow and 
metabolic activity. 

Results 

Gradient Artifact Suppression 

Tested was the principle method of triggered sampling by recording the analog 
25 waveform that drives the gradient coils, to determine the efficacy of the method 

independent of physiological signal fluctuations. The uncorrected signal from a single 
channel appears in Figure 15 (top). Then calculated was the average of 30 repetitions to 
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produce a representation of the gradient activity. The averaging process removes noise 
uncorrected with the imaging gradients. Finally, this averaged signal was subtracted from 
the uncorrected signal, yielding the corrected signal appearing in the lower portion of 
Figure 15. After correction, the large ampEtude gradient activity is removed completely 
5 and all that remains in the signal is the small, uncorrected, noise fluctuation. 

The echo-planar imaging gradient artifact, which lasts for 23 ms and goes through 32 
sinusoidal oscillations in that period, is sampled with only four data points for these data. 
The digitized waveform is thus a very crude record of the actual gradient activity — well 
below the Nyquist frequency. In addition, the morphology of the gradient waveform is 
l o grossly different for each of the nineteen slices that are acquired in this TR period. This 
is due to the feet that the digitization is phase shifted slightly with respect to the readout 
of each slice location. There is also no amplifier "ring down" following the gradient 
pulses, as there is no high pass filter on the amplifier input 

A gitniW experiment was performed, recording electrical potentials from the human 
scalp with a single differential electrode pair. Once the subject was placed into the 
ima ging system, the scanning protocol was similar to that used above: fourteen slice 
locations with a 3 second tr. As acquisition of each slice requires approximately 38 msec 
of gradient activity, the gradients can be expected to obscure about 20% of the EEG 
record. As a consequence of several system non-idealities (such as the small inductance 
of the EEG leads, there will still be some electrical "ring-down." Figure 8 (bottom) 
shows 3 s of the EEG record during scanning, sampled at 200 Hz and triggered by the 
scanner. Prominent gradient artifacts are present even after low pass filt ering . Figure 16 
also shows a series of 20 successive 3 second traces after gradient artifact correction. 

JESQ Energy Mapping 

*» 2 5 The subject was scanned while he lay prone in the magnet with his eyes open. The raw 
EEG data were processed as outlined above, first blanking the gradients and then 
removing the ballistocardiogram. The data were then submitted to a Fourier analysis and 
the energy at each of five pre-defined frequency bands was determined for each TR (that 
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is, fat each image time point). For each frequency band, therefore, it was possible to 
generate a separate time course. 

Shown in Figure 17A is the energy as a function of time in the alpha and theta bands for 
a four and a half minute period while the subject was at rest. The energy levels in these 
5 bands are largely independent for the first minute and a half of this session and seem to 
co-vary for the latter portion of the recording. Figure 17B demonstrates the estimated 
time course of the BOLD signal response based on the (untested) assumption that the 
hemodynamic response related to the EEG is similar to that seen in activation studies 
(e.g., (Cohen 1997; Cohen andDuBois 1999)). 

10 Using the latter as a reference function, then calculated were correlation maps for each 
pixel location with the energy intensity in each previously defined frequency band using 
scanSTAT. The images were first subjected to a simple 9 pixel smoothing to reduce 
pixel noise. The resulting b rain maps are shown in Figures 18A-C, with a separate map 
for the pixel correlation with each of the frequency bands. The images are remarkable in 

15 showing very high correlations, little artifactual activation in white matter or CSF, and 
substantial symmetry, all of which suggest that the artifact content is low. Note also that 
the signal intensity in extra-striate areas shows a substantial signal decrease with increases 
in alpha level 

Conclusions 

2 0 While there have been numerous obstacles to the fusion of functional MEJ and 

electtoencephalogtaphic data, the purely technical challenges of eliminating the cross 
contamination of the recorded signals seem to have yielded to rather sttaightforward 
engineering solutions of somewhat surprising simplicity. The tools outlined here are all 
both easy to construct and made from inexpensive components. The complete software 

25 solution can be made to run as a standalone unit using the LabView pro gramming 

environment The availability of high quality integrated instrumentation amplifiers such 
as the INA114 from Texas Instruments Corporation (Figure 14) makes the analog 
engineering uncomplicated as well. 
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Certain factors affect tie extent to which the data can be interpreted accurately. For 
example, the energetics of the spectral EEG signal are themselves largely unknown. 
Some components, such as the theta rhythm, have long been suspected to be driven by 
subcortical generators and might plausibly be associated with increased thalamic activity 
5 reflected as an increase in blood flow. The alpha rhythm, however, may well be in trinsi c 
to the cerebral cortex (though influenced by thalamus) and is often interpreted as a sort 
of cortical resting state. What is clear, however, is that by continuing to study the 
relationship between scalp (and presumably brain) electrical pot e ntials and magnetic 
resonance signal intensity (presumably indicating increases or decreases in cortical blood 
10 flow) it should be possible to understand better the physiological basis of the 
electroencephalogram (Schomer, Bonmassar etaL 2000). 

The signal processing methods developed to solve the problems in simultaneous EEG 
and fMRI have a broad range of additional applications. For example, these methods can 
be used, without modification, to study electrical evoked potentials (Bonmassar, Anami et 

15 aL 1999) and their localization or, conversely, to yield a better temporal of the "single- 
trial" evoked MR response (Buckner, Bandettini etaL 1996). The clinical applications of 
fully integrated fMRI and EEG are substantial, the most immediate harvest will likely be 
in seizure source localization as an adjunct to other diagnostics used in surgical p lanning 
(Engel 2000). Interictal spikes are a common finding in the epileptic brain yet even when 

2 0 recorded from the cortical surface do not provide definitive and reliable source 

localization. Recendy, however, fast ripples have been reported in rat models and in the 
epileptic human brain, that are associated closely with seizure foci (Bragin, Wilson et aL 
2000). A still speculative interpretation is that such activity represents rapid spontaneous 
action potentials that are ul timat ely propagated and observable as interictal spikes. If so, 

2 5 this gives added hope to the idea that localization by fMRI of interictal spikes might 

become a reliable means to plan resective brain surgery, as has been proposed already by 
several others (e.g., (Warach, Levin et aL 1994; Warach, Ives etaL 1996; Allen, Polizzi et 
aL 1998; Ramabhadran, Frost et aL 1999; Hoffmann, Jager et aL 2000)). 
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Although EEG recorded during MR imaging is particularly noisy, most of the artifacts 
contended with here are present in conventional EEG, though at lower amplitude. For 
example, the activity of the heart, particularly the electrocardiogram, is a contaminant 
that can be corrected in the same manner as the minimization of the ballistocardiogram. 
5 AC line noise is often present as well, and could be eliminated by using a triggered 

sampling approach, timing each block of samples to the power line oscillations — these 
could be detected readily with a phase-locked loop. 

Digital sampling based on precise synchronization with well characterized noise sources 
(Le., triggering on the basis of scanner gradient timing) has applications that extend well 
10 beyond Magnetic Resonance Imaging. It is clear, for example, that a comparable 
approach could be used to remove noise sources associated with AC power line 
oscillations; these include the removal of <4 hum" from digitally processed audio, light 
flicker from fluorescent illumination, and other contaminants of digitally sampled signals. 

Exam ple 3s Tnrr easing Dyna mic T*fln ge Available in Digitized Electro- 
15 encephalography Signals 

The scalp electrical potentials used in the EEG contain both time varying and static (DQ 
components. Often the DC offset is much larger than the EEG, but it is seldom of 
* interest for clinical diagnostic purposes, as it contains essentially no information. 
However, it does cause trouble as it increases the dynamic range needed to digitize the 

2 o EEG signal. For example, the EEG may be only a few microVolts (jiV), while potentials 
of several tens of milliVolts (mV) may exist between electrodes, or as a result of the 
chemical electrode potential developed when the electrode is in contact with the scalp. 
The depth of the signal digitization will be reduced by the ratio of the EEG to the DC 
potential: assuming a DC offset of 10 mV, an EEG signal of 10 jxV, and a 12 bit analog 

25 to digital converter (ADC), the 4096 different levels representable by the ADC will be 
reduced to only 4 levels for the EEG. Clearly this is unacceptable, as the quantization 
noise will dominate the recordings. 
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For this teas on, conventional KEG amplifiers ate equipped with AC-coupled (high pass) 
inputs', usually a capacitor separating the output of the first stage amplifiers from the 
input to the ADC. As the frequencies of interest in EEG can be quite low, the inputs will 
typically have tim e constants of several seconds, allowing fluctuations of 1 Hz or so to 
5 pass without significant attenuation. One consequence of this AC coupling is that it 

creates a tim e constant for signal recovery if the input saturates and, because these filters 
must pass very low frequencies, the settling time for the analog signal to come back to 
the center of its nominal range ca n be quite long. For conventional KEG recordings this 
is acceptable, as large DC shifts are not frequent 

10 In fMRI, the recovery time associated with AC coupling is a particular problem, as the 
gradient-induced artifacts can be large enough to bring the input stages into saturation. 
When this occurs, the amplifiers that follow can be pinned at either the positive or 
negative supply rails for several milliseconds. When the gradients cease, the amplifiers 
may require significant settling time such that the gradient artifact substantially outlasts 

15 the gradient event 

A recent report for example, from Lovblad, et al, reported that even when using a 
reduced gradient activity BURST (Hennig and Hodapp 1993) sequence, "The EEG 
could not be interpreted during the artifact caused by tlie excitation pulse, but the 
recording becomes readable in less than 1 second (approximately 100 msec) after 
2 0 completion of the BURST 11 (Lovblad, Thomas et aL 1999). Thus, if they were to use, for 
example, an eight slice acquisition with a 2 second TR (repetition time), the EEG signal 
would be obscured more than V»a1f of the time by either gradient activity or amplifier 
recovery. 

To remove the DC offset, and to reject saturation problems, the invention provides a 
2 5 circuit that behaves as illustrated in Figure 19. 

In this circuit, the high frequency common mode artifacts typically present in EEG 
recorded in the MRI environment are first attenuated using passive components before 
being differentially amplified. The differential amplification further reduces common 
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mode artifacts. The data ate then presented to a high pass ("anti-aliasing") filter to both 
attenuate the artifacts and to minimise the possibility of aliasing in analog to digital 
conversion. The final stage of this circuit performs the dual role of amplifying the signal 
further and removing any DC bias. This is done by comparing the signal differentially to 
5 a DC reference value sampled from the signal itself. The sampled signal is stored under a 
command here shown diagrammatically as a switch connected to ground. 

Figure 14 shows one preferred implementation, in which the first high pass filter consists 
of a resistor-capacitor (RQ network 20. The first differential amplifier is an integrated 
10 instrumentation amplifier 28. The active high pass filter is a 24 dB Butterworth filter. The 
sample/hold circuit 24 stores the DC signal as charge across a low-leakage capacitor 37 
with a compensation circuit for the non-zero input offset current found in the amplifier 
that drives it. The final differential stage is a standard operational amplifier 34. 

15 With more specific reference to Figure 14, a pair of electrodes is connected to the 

differential input amplifier stage 20, bypassed at high frequencies by parallel capacitance. 
The input amplifier 28 itself may be an integrated circuit, such as an INA114 from Burr- 
Brown corporation. A pair of matched 5 megohm resistors supply bias current to the 
integrated circuit. The output from the amplifier is used as an input to a sample and hold 

2 0 stage 24, which utilizes an integrated circuit 36 such as the UF398 from National 

Semiconductor corporation, which stores the DC voltage across a capacitor 37. A series 
resistance of 200 k£2 provides a time constant of 5 seconds to the sample and hold stage. 
The output of the sample and hold stage is applied to the offset reference pin of the 
differential amplifier component. An offset trim adjustment 40 is provided to adjust for 
25 small static DC offsets in the integrated circuit 36. An input 39 is provided to sample and 
hold device 36. 

The DC corrected output is presented to the input of a low pass filter 22, that utilizes 
operational amplifier integrated circuits 30 and 32, such as the TL072 from Texas 

3 o Instruments corporation. With the component values indicated in Figure 14, the low- 
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pass filter is implemented in a Chebyshev configuration to provide approximately 30 
decibels/octave of attenuation at frequencies above the selected pass band. A final gain 
stage 26 includes a gain trimming adjustment 38 to match the gains across channels. 

5 Example 4: Further Reduction of Setisithrity to Electrical Interference Prom 
Ongoing SgaaslBg Artj Y^tY 

In an additional preferred embodiment, diagrammed in Figure 20, the differential 
recording input stages are supplemented by a shield driver circuit that further reduces the 
sensitivity to electrical interference from the ongoing scanning activity, or other common 

10 mode signals in the leads. In this embodiment, die output from the initial amplifiers is 
coupled through an isolation amplifier to reduce the possibility of circulating currents 
from the RF pulses saturating later components of the amplifiers (Figure 20). The DC 
offset signal is first detected digitally, then held in a digital latch, and returned via a digital 
to analog converter as an offset correction. This offset correction is also coupled 

15 through an isolation amplifier. In this embodiment, die analog output is coupled optically 
through die RF shielded room to minimize any problems of corruption of the MR signal 
through factors such as ground loops. Provision is made for the recording of multiple 
channels of EEG data through the use of multiplexing circuitry. 
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The foregoing description of preferred embodiments of the invention has been 
presented for the purposes of illustration and description. It is not intended to be 
exhaustive or to limit the invention to a precise form disclosed. Many modifications and 
5 variations are possible in light of the above teaching. It is intended that the scope of the 
invention be limited not by this detailed description, but ratifcier by the claims appended 
hereto. 
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WHAT IS CLAIMED IS: 

1 . A method of reducing contamination of electrical signals recorded in the 
presence of repeated interference contamination comprising: 

(a) obtaining an electrical signal, wherein the electrical signal was recorded in 
5 the presence of a contaminating signal; 

(b) detecting a timing signal that occurs at a fixed time point during die 
electrical signal relative to the onset of the contaminating signal; 

(c) digitizing the electrical signal, wherein the digitizing begins with die 
riming signal; 

10 (d) analyzing a plurality of digitized electrical signals, wherein the electrical 

signals are synchronized with respect to the timing signal, to obtain an 
estimated contaminating signal; and 

(e) subtracting the estimated contaminating signal from the digitized 

electrical signal, thereby reducing contamination of the electrical signal. 

15 2. The method of claim 1, wherein the analyzing to obtain the estimate of die 
contaminating signal comprises averaging the elec trical signals. 

3. The method of claim 1, wherein the analyzing to obtain the estimate of the 

contaminating signal comprises calculating a weighted average of the electrical 
signals. 

2 0 4. The method of claim 3, wherein the estimate of the contaminating signal is 
biased towards recent events. 

5. The method of claim 4, wherein the estimate of the con taminating signal 

comprises a first sum consisting of the nth electrical signal added to a scalar 
multiple, w, of a prior estimate of the contaminating signal, divided by a second 
25 sum obtained by adding a series l-fw 2 4-w 3 4-w 4 ... + ... w n . 
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6. The method of claim 1, 2 or 3, wherein the electrical recording comprises an 
electrophysiological signal. 

7. The method of claim 1, wherein the estimate of die contaminating signal is 
multiplied by a scalar prior to the subtracting of step (e). 

5 8. The method of claim 4, wherein the electrophysiological signal comprises an 
electroencephalogtaphic recording, an electtomyelographic record ing , an 
electrocardiographic recording or a measure of galvanic skin resistance. 

9. The method of any one of claims 1-8, wherein the interference comprises 
interference arising from inductively coupled magnetic fields. 

10 10. The method of any one of claims 1 -9, wherein the interference comprises 
interference arising from alternating current (AC) line noise. 

11. The method any one of claims 1-10, wherein the digitizing is performed at a 
sampling rate below the Nyquist rate for the contaminating signal. 

12. The method of any one of claims 1-1 1, wherein the electrical signal obtained in 
15 step (a) is passed through a low pass filter prior to the digitizing, at a frequency of 

approximately one half of the frequency at which the electrical signal is sampled. 

13. The method of claim 12, wherein the low pass filter passes signal frequencies of 
less than about 200 Hz. 

14. The method of any one of claims 1-13, wherein the method is performed 
2 0 concurrently with Magnetic Resonance Imaging of the subject 

15. The method of claim 14, wherein the electrical signal comprises an 
electrophysiological signal and the contaminating signal comprises gradient 
activity. 
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16. The method of claim 15, wherein the electrical signal comprises an 
electrophysiological signal and the con taminating signal comprises radio 
frequency transmitter activity. 

17. The method of any one of claims 1-16, wherein tie digitising is performed at a 
5 rate of about 200 to about 5000 samples per second 

18. The method of claim 1, further comprising removing a DC offset from the 
electrical signal by analog subtraction prior to the digitising. 

19. The method of claim 1 8, wherein the DC offset is measured and subtracted from 
the electrical signal using a differential amplifier. 

10 20. The method of claim 19, wherein the DC offset is measured by analog to digital 
conversion, and averaged over a time period long compared to lowest 
frequencies of interest in the electrical signaL 

21. The method of claim 20, wherein the time period is approximately 10 times 
longer than the lowest frequencies of interest in the electrical signaL 

15 22. The method of claim 20 or 21, wherein the analog subtraction comprises 

converting the averaged signal to an analog voltage and electrically subtracting 
the averaged signal from the electrical signal through differential amplification. 

23. The method of claim 1 9, wherein the DC offset is measured in an analog 
integrator having a time constant long compared with lowest frequencies of 

20 interest in the signaL 

24. The method of claim 23, wherein the time constant is approximately 10 time the 
lowest frequencies of interest 

25. The method of claim 7, wherein the electrophysiological recording comprises an 
electro encephalogram and is recorded concurrently with magnetic resonance 

25 image acquisition. 
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The method of claim 25, wherein the electrophysiological recordings are used to 
inform interpretations of the magnetic resonance imag es. 

The methods of claim 25 or 26, wherein the electrophysiological recordings are 
used in a statistical analysis of change in intensity of the magnetic resonance 
signal. 

The method of any one of claims 25-27, further comprising determining a 
correlation between change in intensity of the magnetic resonance signal and a 
feature of the electrophysiological recording. 

The method of claim 28, wherein the correlation is used to mal ^ statistical 
images that represent an association between the electrical signals and the 
intensity of the magnetic resonance signal intensity. 

The method of claim 28, wherein the feature of the electrophysiological 
recording comprises a time course of signal intensity change in denned frequency 
bands contained in the electrophysiological recording. 

The method of claim 30, wherein the defined frequency bands correspond to 
standard ranges used for clinical interpretations of the electroencephalogram. 

The method of claim 31, wherein the standard ranges are selected from the gtoup 
consisting of from 0 to approximately 4 Hz (the Delta band), from approximately 
4 to approximately 8 Hz (the Theta band), from approximately 8 to 
approximately 12 Hz (the Alpha band), from approximately 12 to approximately 
30 Hz (the Beta band), and from approximately 30 Hz and greater (the Gamma 
band). 

The method of any one of claims 25-32, further comprising convolving the time 
course of the electrophysiological signal with an estimate of the magnetic 
resonance hemodynamic impulse response function. 
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34. A method of reducing magnetic interference during electrophysiological 
recording from a subject by measuring an electrical potential difference between 
a pair of electrodes, wherein the pair of electrodes communicate with a 
differential amplifier via electrical connections, the method comprising twisting 

5 the electrical connections together, thereby reducing magnetic interference* 

35. A method of reducing magnetic interference during electrophysiological 
recording from a subject by measuring an electrical potential difference between 
a pair of adjacent electrodes, wherein each electrode comprises two leads, the 
method comprising twisting each lead together with a lead of an adjacent 

10 electrode, thereby reducing magnetic interference. 

36. The method of claim 34 or 35, wherein the electrophysiological recording 
comprises an electroencephalographic recording. 

37. The method of any one of claims 34-36, wherein the method is performed 
concurrently with Magnetic Resonance Imaging of the subject 

15 38. An apparatus for processing digitized electrical signals in the presence of a 
repeated contaminating signal comprising: 

(a) a signal processor 90 adapted to receive a recording of an electrical signal; 

(b) a detector 92 adapted to detect a timing signal that occurs at a fixed time 
point during an electrical signal relative to the onset of a contaminating 

20 



(c) a signal accumulator 94; and 

(d) a processor 96 adapted to subtract averaged waveforms from an electrical 
signal 
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